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Abstract. 

The equation of the Bohr-Mottelson Hamiltonian with a sextic oscillator potential 
is solved for 7-rigid prolate nuclei. The associated shape phase space is reduced to three 
variables which are exactly separated. The angular equation has the spherical harmonic 
functions as solutions, while the (3 equation is brought to the quasi-exactly solvable 
case of the sextic oscillator potential with a centrifugal barrier. The energies and the 
corresponding wave functions are given in closed form and depend, up to a scaling 
factor, on a single parameter. The O'*' and 2+ states are exactly determined, having 
an important role in the assignment of some ambiguous states for the experimental 
P bands. Due to the special properties of the sextic potential, the model can 
simulate, by varying the free parameter, a shape phase transition from a harmonic 
to an anharmonic prolate /3-soft rotor crossing through a critical point. Numerical 
applications are performed for 39 nuclei: ®®“^°®Ru, 

146-150Nd, 152.154 GcI, 154.156j3y^ 1720s, 180-196pt^ igOjjg ^nd 222^^. The 

best candidates for the critical point are found to be and ^20,i26j^g^ followed 

closely by ^^sxe, ^’’^Os, ^^^Pt and ^^sxd. 
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1 . Introduction 

Even if there are more than 60 years since the Bohr-Mottelson model HH was proposed, 
the interest in its solutions has grown exponentially in the last two decades and this 
happened at least from two main reasons. One is represented by the appearance of its 
cornerstone solutions, called E(5) [3], X(5) [1], Y(5) |5] and Z(5) [ 6 ], which describe 
the nuclei situated in the critical points of the shape phase transitions from spherical 
vibrator to a 7 -unstable system, from spherical vibrator to prolate rotor, from axial 
rotor to triaxial rotor, and from prolate rotor to oblate rotor, respectively. The second 
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reason is due to the increasing number of the experimental data for the quadrupole low- 
lying states and which require interpretation. The multitude of various solutions for 
the Bohr-Mottelson Hamiltonian proposed for this scope were systematically reviewed 

in Refs. [THIS]. 

Soon after the critical point symmetries were proposed, another direction started 
to develop by checking if the 7 -rigid versions of some of these symmetries reveal new 
features for the held of the nuclear shape phase transitions. Therefore, the Z(4) [T3] and 
X(3) [15] solutions were indicated as 7 -rigid counterparts of the Z(5) and X(5) models, 
the hrst corresponding to a triaxial shape and the second to a prolate shape. Z(4) 
is described by the Davydov-Chaban Hamiltonian [16] which emerges from the Bohr- 
Mottelson Hamiltonian by freezing the 7 degree of freedom. Going further and imposing 
axial symmetry, one obtains after quantization in the remaining curvilinear coordinates 
a new Hamiltonian na associated to the X(3) equation. Both Z(4) and X(3) involve an 
inhnite square well potential (ISWP) for the (3 variable, a potential which is well known 
in literature to be a fair approximation for the critical point potential of the spherical to 
deformed shape phase transition. Between the advantages of the ISWP are included the 
simple and the parameter free analytical solution. However, these approaches lack in 
flexibility and consequently in experimental candidates. A question that naturally arises 
is if there are other more flexible potentials which can closely reproduce the properties of 
the ISWP and in the same time present new features, but still lead to analytical solutions 
and depend on very few parameters. The quartic anharmonic oscillator potential 
(QAOP) [IT] and the sextic anharmonic oscillator potential (SAOP) [H] applied for 
X(3) and the quasi-exactly solvable sextic potential [I9] used in connection to Z(4) are 
good examples of this challenge’s importance. 

In the present work the ISWP used in Ref. [I5] for X(3) is replaced by the sextic 
oscillator potential [20]. It is worth to mention here that hrstly the Bohr-Mottelson 
Hamiltonian with sextic potential was solved for 7 -unstable nuclei [ 211122 ], then for 
7 -stable triaxial nuclei [23l[23] and ultimately for 7 -stable prolate nuclei [25]. Other 
numerical applications of these solutions were done in Refs. [26H28] . In the present case, 
the variables are exactly separated and the resulted [3 equation is quasi-exactly solved. 
The energy spectra and the wave functions are given in an analytical form depending, up 
to a scale factor, on a single free parameter. For particular values of the free parameter, 
some terms of the sextic potential as or (3"^ vanish leading to parameter free solutions. 
Due to the special properties of the sextic potential, by varying the free parameter, 
a shape phase transition can be covered from a 7 -rigid prolate harmonic vibrator to 
an anharmonic one, passing through a critical point. The energies and the B{E2) 
transition rates associated to the states 0 + and 2 + of the ground band and of the hrst 
two (3 bands are exactly determined, representing thus reliable reference signatures for 
the identihcation of the possible f3 bands. For the rest of the states an approximation 
is involved whose consistency is validated in respect to the numerical results of the 
previous solutions of the same shape phase space, i.e. X(3), X(3)-/9^ [H], X(3)- f3^ m 
and X(3)-/3® [H]. The solution proposed in the present work is conventionally called 
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X(3)-Sextic, in connection to the X(3) model. Numerical results of X(3)-Sextic, for the 
energy spectra of the ground band and of the hrst two fd bands as well as for several 
E2 transition probabilities, are compared with the available experimental data for the 
nuclei: 100 , 102 ^ 0 , He-iaOXg, 132,134^6, 150,152g^^ 152,154^^, 154,156i3y^ 

^^^Os, i30-i96p|-^ i90]dg and ^^^Ra. A special attention is allocated to the experimental 
evidence of the shape phase transition from 7 -rigid prolate harmonic vibrator to an 
anharmonic one and to the identihcation of the heads of the fd bands, respectively. 

The plan of the present work is the following: X(3)-Sextic solution is described in 
detail in Section II, while in Section III its numerical applications are presented and the 
results are discussed from both theoretical and experimental point of view. Finally, the 
main results of X(3)-Sextic are collected in Section IV. For a better understanding and 
future applications of the X(3)-Sextic solution, an appendix with explicit expressions 
for energies and wave functions is added. 


2. X(3)-Sextic solution 


2.1. Separation of variables 

The Bohr-Mottelson model oia with 7 rigidity (7 = 0) leads to the Davydov-Chaban 
Hamiltonian [16], while in a more particular situation of the 7 -rigid prolate shape 
(7 = 0 °) it takes the following form 


where. 


H = - 


— 


r 1 ^ 02 ^ 1 A ' 

2B [/d^d^^ d/d^ 3 / 32 ^''’'^ 


1 d , 1 

sin 6 89 89 sin^ 9 8(fP' 


V{^). 


( 2 . 1 ) 


( 2 . 2 ) 


The separation of variables is achieved choosing the wave function as T(/3,6',0) = 

F{ld)YL,MA0,4>y 

- Ae,^VpM,(0,0) = L{L + l)ypM,(0,0), (2.3) 


1 d 2 ^ + 1 ) 


7/9 


+ 


+ ^(/9) 


F{/d) = eF{fd), 


(2.4) 


/d^d/d^ did 3^^ 

where by v{/d) = 2BV{(d)/lTf and £ = 2BF/hf are denoted the reduced potential and 
energy. Eq. 02.41) was solved before for an ISWP [TB], a pure harmonic oscillator 
(HO) [ 13 , a QAOP [1^ and a SAOP [18]. In the following, a sextic potential will be 
considered [20] . As can be seen from Table [ 1 ] where a comprehensive account of 7 -rigid 
solutions is presented, the sextic potential contains both fd^ and fd^ terms compared to 
QAOP and SAOP. On the other hand, all these three potentials, for particular values 
of the corresponding free parameters, reproduce closely the ISWP and the HO. 
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Table 1. The potentials in the jd variable and the 7 rigidity values for the relevant 
7 -rigid solutions. 


Solution 

(d potential 


7 

X(3) 

0, if /? < /3a;, 


0° 


cx), if /3 > /3a, 



X(3)-/?2 

~/32 


0° 

X{3)-/3^ 

~/3^ 


0° 

X(3)-/36 

~/36 


0° 

QAOP 

|q:i/3^ + 0 : 2 / 3 "^, 


0° 


ai > 0, (72 > 0 



SAOP 

^(71/3^ + 0 : 2 / 3 ®, 


0° 


0(1 > 0, 0(2 > 0 



X(3)-Sextic 

(62 - 4ac)/32 + 206/34 

+ o2/3®. 

0° 


c, a > 0, b E R 



Z(4) 

0, if /3 < /3a;, 


CO 

0 

0 


00 , if /3 > /3a, 



Z(4)-/?2 

~/32 


30° 

Z(4)-Sextic 

(62 - 4ac)/32 + 206/34 

+ o2/3®. 

30° 


c, 0 > 0, b E R 




2.2. Solution of the j3 equation with sextic potential 


In what follows, it is preferable to write Eq. fl2.4p in a Schrodinger form. This is realized 
by changing the wave function as F{/3) = 


L{L + 1) 
d ^^ 3/32 


+ v{/3) 


(p{/3) = 


(2.5) 


Further, Eq. (12.bh is brought to a form associated with the quasi-exactly case of the 
sextic oscillator potential with a centrifugal barrier [ 20 ] by considering the relations: 


L{L + 1) 
3 



( 2 . 6 ) 


7 / 3 ) 



(. + 1 + m) 


+ 2abfd^ + M 


0,1,2,... (2.7) 


Here by quasi-exactly solvable one understands that the eigenvalue problem is exactly 
solved only for a part of the spectrum, namely for M+1 eigenstates. The sextic potential 
fl2.7p depends on two parameters, a and b, and on an integer M, but also on the quantum 
number L through s. Apparently, the potential is state dependent, but, as it will be 
shown, M can be used instead to keep the potential independent on L by forcing that: 

s + - + M = const. = c. 


( 2 . 8 ) 
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s{L) 



L{L + 1) 



(2.9) 


it is observed that the condition (12.81) is impossible to be fulhlled for L > 4, s taking 
irrational values while M is an integer number. Instead, for L = 0 and 2 the expression 
of s given by Eq. fl2.9p reduces to a simpler form: 


,'(L) = (2.10) 

By comparing, in Table [2l s(L) and s'{L) for different values of L it can be appreciated 
that the latter represents a good approximation of the former at least for L < 10. Using 


Table 2. A comparison of s{L) and s'{L), given by Eqs. (12.911 and (12.1011 for even 
L < 10. 


L 

0 

2 

4 

6 

8 

10 

s 

0.75 

1.25 

1.81 

2.39 

2.96 

3.54 

s' 

0.75 

1.25 

1.75 

2.25 

2.75 

3.25 


s'{L) even for L > 4, the condition (12.8p becomes 

L 5 
c = M + - + -, 

4 4 

which is constant if M decreases by one unit while L increases by four units: 

5 
4 


(M, L) : (A', 0), (A - 1,4), (A - 2,8 ),^ A' + ^ = cf\ 


( 2 . 11 ) 


( 2 . 12 ) 


(M, L) : (A, 2), (A - 1,6), (A - 2,10),... ^ A + 1 = cf >. 


(2.13) 

Here by K is denoted the maximum value of M. Finally remain two potentials 
independent on L, one for L = 0,4,8,... and other for L = 2,6,10,..., which slightly 
differ through the coefficient of by the quantity = 1/2. The sextic potential 

equation (12.5p can be more simplihed by reducing the number of parameters through 
the change of variable (3 = and adopting the notations a = h/y/a and Sy = ejyfa\ 

L{L + 1) 


,{K) 


dy^ 


3y^ 


(y) 


viy) = ^yviy), 


(2.14) 


where. 


vl^\y) = - 4cif^)i/^ + 2ay^ + / + m = 0, 2. (2.15) 

(K) 

Here, Um ^ are constants added in order to minimize the slight difference between the 
two potentials fl2.15p . Thus, Um^ are hxed such that the potentials to have the same 
minimum energy: 
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Figure 1. The shape evolution of the scaled sextic potential v^\y), given by Eq. 

(I2.15|l . as a function of the free parameter a for = 13/4 and = 0. 

where by ?/q'^ are denoted the extreme points of the scaled potential (12.151) : 

a2 + 12c^^J , m = 0,2. (2. 

An elegant representation of the shape evolution of the sextic potential fl2.15p as a 
function of the free parameter a is shown in Fig. [H Thus, coming from a ^ oo towards 
0 = 0 through the point Oc, a shape phase transition is covered from a harmonic vibrator 
to an a n h armonic one crossing a critical point for which the potential is flat. Going 
farther to a ^ —oo a simultaneous spherical and deformed minima appear separated 
by a maximum. Actually this representation is a projection of the sextic potential 
evolution from the plane of parameters (a, b) [22] on the real axis of the free parameter 
a. Indeed, in Ref. [22] the critical point corresponds to the parabola a = b'^/Ac, while 
here is given by = 4c, where a = b/y/a. Another important remark is that points 
ac = ±2y^ are moving away from origin as the K value increases. 

So far, all the above efforts were done to bring Eq. (12.bh with a sextic potential [20] 
to a solvable equation for the y-rigid prolate nuclei and then to have a more convenient 
form of it by reducing the number of parameters. Further, the attention will be focused 



( 0 " = 5 (- 2 “±'/ 

























Sextic potential for ■y-rigid prolate nuclei 

on the solution of Eq. fl2.14p . Therefore, choosing as an ansatz the function 


V 




2 g 4 2 


( 2 . 18 ) 


Eq. (12.1411 with the potential fl2.15p is reduced to the equation 


/ ^ 4s' — Id 


y dy 


dy 


d 

dy 


2M 


P(^\y^) = 


(2.19) 


Vdj/2 

= APW(l/2) 

where P^^l are polynomials in y‘^ of order M. Writing Eq. fl2.19p in a matrix form as 
in the Appendix of Ref. [23] , the eigenvalues A are obtained as 

A = yW = £ _ 4as' - _ (— 

where the last term, containing the mean value of ?/^, was extracted from the centrifugal 
term of Eq. fl2.14p such that s(L) could be approximated by s'(L) for L > 4. By adding 
this term to the energies (12.2011 . arises the possibility to describe very well also the states 
with > lO’*', for which s(L) and s'{L) otherwise start to diverge. Finally, the total 
energy of the system is easily deduced from Eq. (I2.20p : 


( 2 . 20 ) 


En,L — 


2B 


+ 3) + u^^\a) + 


1 L 

( 2 / 2 ) 6 


1-1 

2 


( 2 . 21 ) 


The total energy (I2.2ip depends on two quantum numbers, the order of the wave 
function’s zero n and the total angular momentum L, on a hxed integer K, on a scale 
factor {h^^fa) / {2B) and on a free parameter a. There are two possibilities for numerical 
applications to the experimental data. The hrst one is to normalize the energies (I2.2ip 
to the ground state energy and to £t both the scale factor and the free parameter or 
the second one by normalizing to the hrst excited state energy, as 

En,L - Po,o j - 2 . 22 ) 


En,L — 


Po,2 — Po,0 


and then to £t only the free parameter. In the present paper, the second method is 
adopted with n = 0, n = 1 and n = 2 corresponding to the ground band, the hrst and 
second (3 bands, respectively. 


2.3. E2 transition probabilities 

Another important aspect of the present approach comparing to the previous ones, i.e. 
QAOP [T7| and (SAOP) [H] is that its associated quasi-exactly solvable method allows 
to determine the expressions for the wave functions. In this way, numerical results for 
B{E2) transition rates can be performed. The total wave function is a product of the 
radial and angular components: 


(2.23) 
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where YL^ML{d,4>) are the well known spherical harmonics, while are expressed 

in terms of the normalized ansatz functions (12.181) . The harmonic transition operator 
for 7 -rigid prolate case gets the simplihed form [T5] : 




(2.24) 


where t is a scale parameter which is dropped for normalized B {E2) values. Therefore, 
the dual expression of the reduced E2 transition probabilities, normalized to the 
transition from the hrst excited state to the ground state, is: 


T 

-L 1 


n,L,n',L' — 


B{E2] n,L ^ n', L') 

fi(h;2;0,2 ^ 0,0) 


(~iL2L' j 
^ 000 ^n,L-,n',L' 
r^220 j 
O' 000-'0,2;0,0 


(2.25) 


where the radial matrix element In,L-,n',L' can be given either in /? or y variable: 


Ln,L-,n',L' — 




(M) 




(M) 


(2.26) 


= a-'/" 


v^^L{y)yVnyL'{y)^y- 


(M) 


As it will be shown further, X(3)-Sextic presents unique properties from both 
theoretical and experimental point of view by comparison with the previous similar 
solutions, making it a powerful tool in the description of the 7 -rigid prolate nuclei. 


3. Numerical results and discussion 
3.1. Theoretical aspects of X(3)-Sextic 

In the following, the entire discussion will be restricted to the case K = 2 for which, 
as can be deduced from Eq. fl2.13p . the properties of the most representative low-lying 
states of the ground band (up to lO'*'), the hrst ft band (up to 6 ''') and the second (3 band 
(up to 2 +) are reproduced. At this point, the index K will be dropped as it has a hxed 
value. If necessary, results for K > 2 can be considered without any problems [23]. The 
beauty of the case A' = 2 is that explicit analytical expressions for energies and wave 
functions are still possible. For a better understanding and possible future applications 
of X(3)-Sextic solution, the cases corresponding to M = 0, M = 1 and M = 2 are 
detailed in the Appendix. A plot of the energies given by Eq. fl2.22p for K = 2 and 
a G [—10,10] are presented in Fig. [2] from which one can draw some useful insights. 
For a —)■ 00 , X(3)-Sextic reproduces the results of the 7 -rigid prolate harmonic vibrator 
called X(3)-/9^ [T7], while for a close to zero a 7 -rigid prolate anharmonic vibrator 
shows up. Special situations appear for a —)■ — 00 , where the ground band energies 
maintain hnite values while the (3 band energies go to inhnity, and for a G [2^/cq, 2 yT 2 ] 
(the narrow vertical area) where a critical point of a hrst order shape phase transition 
occurs. In the region a G [—2yTo, — 00 ], as can be seen from Fig. [H simultaneous 
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spherical and deformed minima appear separated by a high barrier. On the other hand, 
in the interval a G [2y/co, 2 y^] the potential is completely flat as an inhnite square 
well potential, which is known from the literature that usually corresponds to critical 
points. Indeed, by analyzing the hrst derivative of the energy fl2.22p in respect to 
the free parameter a it is found that it has a discontinuity for an ac G [ 2 y/co, 2 y/c 2 \. 
The interpretation of this result is that a hrst order shape phase transition takes place 
between a 7 -rigid prolate harmonic vibrator and a 7 -rigid prolate anharmonic vibrator. 
Actually this is not a surprise, similar results for different shape phase transitions being 
reported in Refs. [8l[T9l|22l|3Tl[32]. An analogy with X(5), which describes the critical 
point of a hrst order shape phase transition from spherical vibrator to 7 -soft prolate 
nuclei, can be also made. 

By analyzing the expression of the sextic potential (I2.15p it can be observed that 
the terms y'^ and y‘^ are cancelling for a = ±2{cm^y^‘^ and a = 0 , respectively, leading 
to free parameter solutions. The energy spectra fl2.22p and some E2 transitions in these 
particular cases are presented in Fig. [3] for K = 2. If the energy spectrum for a = 0 
doesn’t show nothing out of common, that for ac is quite interesting from the point 
of view of its structure. Firstly, the states are grouped two by two in each band, and 
secondly an approximate degeneracy appears for states of diherent angular momenta 
(AL = 4) belonging to diherent bands. Also the B {E2) transitions are more intense for 
a = ac than a = 0, while the energies are higher for a = 0 than those for a = ac- An 
interpretation of the results at ac in terms of some symmetries would be an interesting 
subject, especially that the dynamical symmetry of X(5) is still unknown. It should be 
noted that degeneracies in the critical regions have been also indicated in Refs. [331134] . 
while an Euclidian dynamical symmetry unifying several of them has been proposed in 
Ref. [35]. 

In order to verify the consistency of the approximation fl2.10p . several energy ratios 
(I 2 . 22 P are compared in Table [3] with those corresponding to X(3)-/3^ [IT], X(3)-/d^ [T7] . 
X(3)-/3® [18] and X(3) [15], respectively. The resulted root mean square (rms) value 
is increasing as one goes to higher anharmonic solution. As was expected, the results 
of X(3)-Sextic are closer to those of X(3)-/l^, X(3)-/5"^, and X(3)-/9® because it contains 
all these terms, while for a hdel reproduction of the ISWP the inclusion of higher 
anharmonicities is necessary [36] . 


3.2. Applications of X(3)-Sextic to experimental data 


The X(3)-Sextic solution is applied to describe the available experimental data for the 
ground and the hrst two jd bands of 39 nuclei, namely, ^^®“^^°Xe, 

132,134Ce, 146-150Nd, 150,152g^^ 152,154Qd, 154,156i3y^ 1720s^ 180-196p^^ igOpg 222Ra. 


The comparison for the energy spectra is given in Tables ID and [5] while for the E2 
transitions is shown in Tables [ 6 ] and d All these data are htted involving only the 


free parameter a, its values and the corresponding rms being indicated in the same 
tables. The agreement of X(3)-Sextic with the experimental data is good in general 
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Figure 2. The energy spectra of the ground band and of the first two /3 bands, 
given by Eq. (12.221) . are plotted as a function of the free parameter a in the interval 
a S [—10,10] for K = 2. The energy lines are indexed by where n = 0 for the 

ground band, n = 1 for the first /3 band and n = 2 for the second /3 band. 


Table 3. Several energy ratios (I2.22|) given by X(3)-Sextic are compared with the 
predictions of the parameter free models X(3)-/3^ [17], X(3)-/3^ [17], X(3)-/3® [18] and 

x(3) m- 


Model 

R-0,4 

R-0,6 

R-0,8 

O 

o 

R.1,0 

Ri,2 

R.1,4 

Ri,6 

R2,0 

R2,2 

rms 

a 

X(3)-d" 

2.13 

2.11 

3.27 

3.24 

4.42 

4.40 

5.58 

5.56 

2.00 

2.02 

3.00 

3.02 

4.13 

4.11 

5.27 

5.21 

4.00 

4.05 

5.00 

5.06 

0.04 

25.6 

X(3)-d‘‘ 

2.29 

2.45 

3.72 

3.72 

5.28 

5.49 

6.94 

6.89 

2.37 

2.48 

3.70 

3.67 

5.29 

5.36 

6.98 

6.76 

5.28 

5.40 

6.83 

6.73 

0.13 

4.92 

X(3)-/3® 

2.34 

2.65 

3.91 

3.99 

5.65 

6.08 

7.57 

7.61 

2.56 

2.76 

4.08 

4.05 

5.94 

6.05 

7.97 

7.59 

6.01 

6.14 

7.93 

7.63 

0.24 

4.04 

X(3) 

2.44 

2.62 

4.23 

4.29 

6.35 

6.81 

8.78 

8.87 

2.87 

3.40 

4.83 

5.19 

7.37 

7.52 

10.29 

9.72 

7.65 

7.97 

10.56 

10.25 

0.35 

1.92 


for most of the considered nuclei and in particular for ^^^Ce, 

^^®Nd, isogg for which small rms values are obtained. Instead, the 

experimental data of ^°®Ru, ^^®Xe, ^^^Ce, ^^°Nd, ^^^Os and ^®^Pt 

are well reproduced only partially, either the band heads or only the (3 bands being well 
described. For the states where the experimental data are of uncertain band, a selection 
has been done based on the corresponding X(3)-Sextic results. This is mostly done for 
the band head state of the second jS band and then for states of the hrst (3 band. These 












Sextic potential for 'y-rigid prolate nuclei 


11 


(2,2) 16.596 


(10,0) 10.327 


2.123 


( 8 , 0 ) 


7.371 


2.012 

(6, 1) I 4.526 



1.761 / = 

S c 

i c 


(4, 1) y 2.343r 




1.557 




(2, 2) y 1 


' 

' 

(0, 2) i 0 

' 

' 


g(n = 0) 


P\ {n = 1) 
a = 0 


h (" = 2) 


(10,0) 8.010 (6, 1) 8.080 (2,2) 8.165 


4.324 

(8, 0) i 6.401 

I 

3.766 
(6, 1) \ 4. 


2.996 


(4, 1) I 6.442> (^ 


2.424 

(2, 2) t 4.2 73 


(2, 2) i 


(0,2) y 0 

g(n = 0) 


O (N 


j6l (n = 1) 


( 0 , 2 ) 


1.323 
6.569 


r- 

00 

00 


04 


O 

1 . 

1 


(6, 1) t 4.140 .If’ ▼ 4.273 _ 

I Il.i52 I 

(4,l)’t 2.745/' t 2.916 _ 

-rt \D 

\ r- m 

2.180 z; 


Pi (n = 2) 


a = 0;^= 3.647 


Figure 3. Energy spectra (12.221) and B{E2) transitions (12.251) for a = 0 (left side) 
and a = ac (right side). The states are indexed by {L,K). 


results can be very useful to clarify the nature of these states. Also, an experimental 
energy ratio of 2.76 is attributed to the band head state of the hrst ft band of ^°^Ru, 
for which X(3)-Sextic predicts the value 2.92. As a matter of fact, the ^°^Ru nucleus 
is important here because its corresponding a = 3.62 value is very close to the critical 
point ttc = 3.65. As has been discussed in the previous subsection, this value of a 
leads to a parameter free description. In the same category are also included the nuclei 
i 2 o,i 26 ^g and ^^^Ra, the last one being close to the a = 0 case. Another interesting study 
carried out here has as starting point the fact that X(3)-Sextic can simulate the y-rigid 
version of X(5). Therefore the experimental data of the most representative candidate 
nuclei for X(5), namely, ^^^Sm, ^^^Gd and ^^®Dy are analyzed in the frame of 

the present solution. It is found that the ground band of these nuclei prefers a y-stable 
structure, while the situation of the hrst two (3 bands is changed in a clear favor for 
y-rigidity. This behavior, also reported in connection to X(3) [I5], rehects the fact that 
a combined scheme of y-stable and y-rigid would be more appropriate for these nuclei. 
Actually, an important step toward this direction has already been done in Ref. EH. 
where a free parameter manages the y-rigidity. 
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Concerning the B{E2) transition rates, given in Tables E] and [TJ the agreement 
of X(3)-Sextic with the corresponding experimental data is good taking into account 
the fact that the free parameter a was htted only for the energy spectra. As can be 
seen, there are some discrepancies in the ground band of some nuclei for which X(3)- 
Sextic predicts increasing values in respect to L while experimental data manifest a 
decreasing trend. This problem can be partly removed by including anharmonicities 
in the transition operator fl2.24p as in Ref. |25]. With this, the B{E2) results between 
states from different bands can also be improved. 

Another important aspect that can be investigated by X(3)-Sextic is whether or 
not a shape phase transition takes place within an isotopic chain. For this purpose, 
the htted values of a are plotted in Fig. [Has a function of the neutron number N for 
the most numerous isotopic chains considered in the present paper, namely Ru, Xe, 
Nd and Ft. In Fig. IH the regions above and bellow the dashed line indicating the 
critical point are equivalent in Fig. |2]with the right side and the left side, respectively, 
of the critical region represented by the gray narrow area. The isotopes of Ru and Nd 
have a similar behavior crossing once the dashed line indicating the critical point and 
having the lightest isotopes and the heaviest ones situated above and below the line, 
respectively. By contrary, the isotopes of Xe cross twice the critical point, once from the 
lightest isotopes towards the medium ones ( 54 ®Xe 62 —>■ s^^Xeeg) and second time from 
the heaviest isotopes towards the medium ones ( 54 °Xe 76 —)■ g^^Xegs)- If the points {a, N) 
of the Xe isotopes are interpolated by a polynomial function, then the hrst derivative 
of the resulted function a = f{N) would present a discontinuity at the point of g^^Xegs 
where the two separated shape phase transitions converge. This last observation can be 
interpreted as follows. If g^^Xegg and g 4 ®Xe 72 are candidates for the critical points of the 
two shape phase transitions ( 54 ®Xeg 2 —)■ g^^Xegg) and ( 54 °Xe 7 g —)■ g^^Xegg), respectively, 
the isotope g 4 ^Xegg plays the role of a critical point of a transition between the two 
transition arms of the Xe isotopes. In other words the transition from gf Xegg to 54 *^X 070 
or g 4 ®Xe 72 is not smooth at all but rather is like a jump. This is not the case for the Ru 
isotopes or for those of Nd where the transition is smooth and unidirectional. Instead, 
a similar jump is observed for the isotopes of Ft, namely between ^I'^Ftiog and yg^Ftiog. 
Nevertheless, because no crossing of the critical point takes place for these isotopes, 
the jump is associated to a critical point of a local shape phase transition taking place 
between the group of more deformed isotopes 7 g°~^^^Ftio 2 -io 6 and the less deformed 
group of 78 “ Ftio 8 -ii 8 nuclei. Concerning the best candidates for the critical point 
of the 7 -rigid prolate harmonic vibrator to 7 -rigid anharmonic vibrator shape phase 
transition, by far it is the 44 ^Rugo nucleus followed closely by g 4 °’^^®Xeggj 2 and gg^Ndgs- 
Surprisingly, 44 "^Rugo was found to be a good candidate also for the critical point of the 
spherical to 7 -unstable shape phase transition [22l[7T] and for a situation intermediate 
between the triaxial and the 7 -unstable limit [72II74] . but what is interesting in the frame 
of the present description is that the isotope 44 '^Rugo falls exactly in the critical point 
of the new transition. Other good candidates for the critical point may be considered 
54 ®Xe 74 and even ^g^Ftng and ^^Osgg despite the fact that no crossing of the critical 
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Table 4. The energy spectra for the ground band and the first two /3 bands given 
by Eq. (I2.22p are compared with the available experimental data [551 - 155] of the nuclei 
98-108ru^ 100,102^0, ii^-isoxe and 132.134^3^ gj.gt Qf nucleus are 

given the experimental data, while in the second line are the corresponding theoretical 
results. In brackets are indicated possible candidate energy states for the corresponding 
predicted data, which were not included in the fit. 


Nucleus 

Ro,4 

Ro,6 

Ro,8 

Ro,io 

Ri,o 

Ri,2 

Ri,4 

Ri,6 

R2,0 

R2,2 

rms 

a 

44 IIU 54 

2.14 

2.21 

3.41 

3.38 

4.79 

4.73 

5.97 

2.03 

2.15 

3.21 

4.49 

5.68 

4.46 

5.57 

0.08 

8.44 

i^°Ru56 

2.27 

2.58 

3.85 

3.90 

5.67 

5.88 

7.85 

7.36 

2.10 

2.66 

3.92 

5.82 

7.31 

5.89 

7.32 

0.37 

4.28 

ifRuss 

2.33 

2.50 

3.94 

3.79 

5.70 

5.64 

7.23 

7.08 

1.99 

2.55 

3.77 

5.54 

6.98 

5.59 

6.96 

0.28 

4.63 

44'^Ru60 

2.48 

2.74 

4.35 

4.14 

6.48 

6.41 

8.69 

8.02 

(2.76) 

2.92 

4.23 

4.28 

5.81 

6.45 

8.10 

6.58 

8.19 

0.40 

3.62 

ifRu62 

2.66 

2.55 

4.80 

4.43 

7.31 

7.14 

10.02 

9.55 

3.67 

4.10 

6.20 

8.55 

11.19 

9.36 

12.20 

0.34 

1.16 

ifRu64 

2.75 

2.49 

5.12 

4.50 

8.02 

7.30 

11.31 

9.92 

4.03 

4.75 

7.01 

9.31 

12.22 

10.43 

13.65 

0.83 

0.73 

1§°M058 

2.12 

2.20 

3.45 

3.37 

4.91 

4.72 

6.29 

5.95 

1.30 

2.14 

2.73 

3.20 

4.47 

5.66 

4.44 

5.54 

0.43 

8.63 

42^M060 

2.51 

2.67 

4.48 

4.19 

6.81 

6.57 

9.41 

8.39 

2.35 

3.09 

3.86 

4.65 

6.90 

8.80 

7.16 

9.07 

0.63 

2.65 

^fXe62 

2.33 

2.28 

3.90 

3.48 

5.62 

4.97 

7.45 

6.26 

2.25 

2.58 

3.35 

3.96 

4.76 

5.38 

6.02 

4.75 

5.93 

0.72 

6.62 

^fXe64 

2.40 

2.61 

4.14 

3.94 

6.15 

5.98 

8.35 

7.48 

2.46 

2.71 

3.64 

3.98 

5.13 

5.93 

7.44 

(5.10) 

6.01 

7.47 

0.49 

4.16 

120v„ 

54 2^666 

2.47 

2.73 

4.33 

4.11 

6.51 

6.32 

8.90 

7.91 

2.82 

2.88 

3.95 

4.21 

5.31 

6.34 

7.94 

(6.93) 

6.45 

8.01 

0.57 

3.79 

^fXees 

2.50 

2.63 

4.43 

4.26 

6.69 

6.74 

9.18 

8.73 

3.47 

3.29 

4.51 

5.02 

7.33 

9.44 

(7.63) 

7.72 

9.89 

0.30 

2.11 

M^Xeyo 

2.48 

2.70 

4.37 

4.17 

6.58 

6.50 

8.96 

8.23 

3.58 

3.01 

4.60 

4.48 

5.69 

6.70 

8.50 

(6.70) 

6.91 

8.69 

0.53 

2.99 

54®Xe72 

2.42 

2.72 

4.21 

4.09 

6.27 

6.29 

8.64 

7.86 

3.38 

2.86 

4.32 

4.18 

5.25 

6.30 

7.89 

(6.57) 

6.40 

7.95 

0.55 

3.83 

^fXe74 

2.33 

2.66 

3.92 

4.01 

5.67 

6.12 

7.60 

7.65 

3.57 

2.78 

4.52 

4.07 

6.10 

7.65 

(5.87) 

6.19 

7.69 

0.44 

4.00 

54°Xe76 

2.25 

2.41 

3.63 

3.67 

5.03 

5.37 

6.75 

3.35 

2.43 

(4.01) 

3.60 

(4.53) 

5.23 

6.59 

5.26 

6.55 

0.50 

5.18 

ii"Ce74 

2.64 

2.65 

4.74 

4.23 

7.16 

6.65 

9.71 

8.56 

3.56 

3.18 

4.60 

4.83 

5.94 

7.11 

9.11 

7.43 

9.47 

0.70 

2.36 

^rce76 

2.56 

2.61 

4.55 

4.30 

6.87 

6.82 

9.09 

8.89 

3.75 

3.42 

4.80 

5.23 

7.56 

9.77 

8.01 

10.31 

0.26 

1.89 
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Table 5. The same as in Table |4j but for the available experimental data [54117(1] 
of the nuclei 150,152gj^^ 152,154 qj^ 154,156i3y^ 172Qg^ 180-196p^^ 190jjg 

222Ra. 


Nucleus 

R-0,4 

Ro.6 

R-o.s 

R.0,10 

R-1,0 

Ri.2 

Ri,4 

R-1,6 

R-2.0 

R2,2 

rms 

a 


2.30 

2.34 

3.92 

3.56 

5.72 

5.14 

7.32 

6.46 

2.02 

2.32 

2.87 

3.45 

3.85 

4.96 

6.26 

4.96 

6.19 

0.64 

5.88 

If^dss 

2.49 

2.61 

4.24 

3.95 

6.15 

5.98 

8.19 

7.49 

3.04 

2.71 

3.88 

3.98 

5.32 

5.94 

7.12 

7.46 

(5.30) 

6.02 

7.48 

0.39 

4.15 

^E°Nd90 

2.93 

2.49 

5.53 

4.50 

8.68 

7.30 

12.28 

9.92 

5.19 

4.75 

6.53 

7.01 

8.74 

9.31 

11.83 

12.22 

(13.35) 

10.43 

13.65 

1.10 

0.73 

M°Sm88 

2.32 

2.36 

3.83 

3.60 

5.50 

5.23 

7.29 

6.57 

2.22 

2.36 

3.13 

3.51 

4.34 

5.06 

6.31 

6.39 

(3.76) 

5.07 

6.33 

0.41 

5.58 

ei^Smgo 

3.01 

2.55 

5.80 

4.43 

9.24 

7.14 

13.21 

9.55 

5.62 

4.10 

6.65 

6.20 

8.40 

8.55 

10.76 

11.19 

8.89 

9.36 

10.62 

12.20 

1.59 

1.16 

^fGd88 

2.19 

2.10 

3.57 

3.22 

5.07 

4.37 

6.68 

5.52 

1.79 

2.00 

2.70 

3.00 

3.72 

4.06 

4.85 

5.16 

3.04 

4.00 

3.83 

5.00 

0.68 

100 

^fGdgo 

3.01 

2.53 

5.83 

4.46 

9.30 

7.21 

13.3 

9.70 

5.53 

4.33 

6.63 

6.50 

8.51 

8.84 

11.10 

11.58 

9.60 

9.75 

11.52 

12.74 

1.51 

0.99 

66 ^Dy88 

2.23 

2.10 

3.66 

3.22 

5.23 

4.37 

6.89 

5.52 

1.98 

2.00 

2.71 

3.00 

3.74 

4.06 

4.96 

5.16 

3.16 

4.00 

4.16 

5.00 

0.67 

100 

efDygo 

2.93 

2.53 

5.59 

4.46 

8.82 

7.20 

12.52 

9.68 

4.90 

4.31 

6.01 

6.46 

7.90 

8.81 

10.43 

11.53 

10.00 

9.71 

12.68 

1.28 

1.01 

re^Osge 

2.66 

2.65 

4.63 

4.00 

6.70 

6.10 

8.89 

7.63 

3.33 

2.77 

3.56 

4.06 

5.00 

6.07 

6.81 

7.62 

6.16 

7.66 

0.77 

4.02 

180p^ 

78 ^^102 

2.68 

2.56 

4.94 

4.40 

7.71 

7.07 

10.93 

9.41 

3.12 

3.92 

5.62 

5.96 

8.15 

8.32 

10.77 

10.86 

(7.69) 

9.03 

11.75 

0.69 

1.31 

^rPtl04 

2.71 

2.57 

5.00 

4.40 

7.78 

7.05 

10.96 

9.37 

3.22 

3.88 

5.53 

5.90 

8.00 

8.26 

10.64 

10.77 

(7.43) 

8.95 

11.64 

0.72 

1.35 

184p^ 

78 ^^106 

2.67 

2.57 

4.90 

4.38 

7.55 

7.01 

10.47 

9.29 

3.02 

3.79 

5.18 

5.77 

7.57 

8.13 

11.04 

10.59 

8.77 

11.39 

0.65 

1.44 

186pf 

78 ^tio8 

2.56 

2.68 

4.58 

4.18 

7.01 

6.53 

9.70 

8.30 

2.46 

3.04 

4.17 

4.56 

6.38 

6.80 

8.36 

8.64 

7.03 

8.87 

0.62 

2.82 

^fPtllO 

2.53 

2.65 

4.46 

4.22 

6.71 

6.64 

9.18 

8.53 

3.01 

3.17 

4.20 

4.80 

7.08 

9.06 

7.39 

9.41 

0.38 

2.40 

190pf 

78 ^tii2 

2.49 

2.69 

4.35 

4.17 

6.47 

6.50 

8.57 

8.24 

3.11 

3.01 

4.07 

4.50 

6.72 

8.52 

(5.65) 

6.93 

8.72 

0.25 

2.96 

^fPtll4 

2.48 

2.65 

4.31 

4.23 

6.38 

6.66 

8.62 

8.57 

3.78 

3.19 

4.55 

4.84 

7.13 

9.14 

7.45 

9.50 

0.30 

2.34 

194px 

78 ^^116 

2.47 

2.64 

4.30 

4.24 

6.39 

6.69 

8.67 

8.63 

3.23 

3.86 

4.60 

4.91 

7.20 

9.25 

7.55 

9.65 

0.32 

2.25 

196pf 

78 ^^118 

2.47 

2.72 

4.29 

4.15 

6.33 

6.45 

8.56 

8.13 

3.19 

2.96 

3.83 

4.39 

6.59 

8.31 

6.76 

8.46 

0.33 

3.25 


2.50 

2.36 

4.26 

3.59 

5.21 

6.56 

3.07 

2.36 

3.77 

3.50 

4.74 

5.05 

6.03 

6.37 

5.06 

6.31 

0.46 

5.62 

irRai34 

2.71 

2.25 

4.95 

4.48 

7.59 

7.29 

10.56 

10.37 

8.23 

7.44 

9.22 

9.90 

11.63 

15.14 

14.19 

18.33 

0.52 

-0.38 
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Table 6. Several B{E2) transition rates normalized as in Eq. (j2.25l) are compared with 
the available experimental data [38115^ of the nuclei 

I32,i34(^g^ In the first line of each nucleus are given the experimental data, while in 
each second line are the corresponding theoretical results. 


Nucleus 

To, 4 , 0,2 

To.6,0,4 

To,8, 0,6 

To, 10 . 0,8 

7 "!, 0 , 0,2 

7 " i , 2 , 0,2 

7 ’ i , 2 , 0,4 

7 ’ i , 2 , 1,0 

44 R-U 54 

0.38(11) 

2.34 

0.40(8) 

3.48 

0.08(2) 

4.62 

5.62 

3.21 

0.21 

1.47 

1.52 

IfRuse 

1.43(11) 

2.23 

4.78(5) 

3.08 

3.98 

4.63 

0.98(14) 

2.83 

0.18 

1.02 

1.24 

lfRu58 

1.48(25) 

2.25 

1.52(56) 

3.15 

1.26(43) 

4.08 

1.28(47) 

4.77 

0.78(14) 

2.90 

0.19 

1.09 

1.28 

44‘*Ru60 

1.43(16) 

2.18 

2.94 

3.76 

0.43(5) 

4.31 

2.65 

0.17 

0.88 

1.15 

44®Ru62 

1.75 

2.07 

2.47 

2.67 

1.22 

0.09 

0.24 

0.69 

44®Ru64 

1.65(20) 

1.66 

1.93 

2.27 

2.43 

0.95 

0.08 

0.16 

0.64 

1§°Mo58 

1 . 86 ( 11 ) 

2.35 

2.54(38) 

3.49 

3.32(49) 

4.64 

5.64 

2.49(12) 

3.21 

Ri 0 

0.21 

0.97(49) 

1.48 

0.38(11) 

1.53 

42^Mo60 

1.20(28) 

2.05 

2.64 

3.31 

0.95(42) 

3.71 

2.23 

0.15 

0.62 

0.97 

^i®Xe62 

1.75(11) 

2.32 

1.58(15) 

3.38 

4.45 

1.56(29) 

5.34 

3.12 

0.20 

1.35 

1.44 

ii®Xe64 

1 . 11 ( 6 ) 

2.22 

0.88(23) 

3.06 

0.49(18) 

3.95 

>0.7 

4.58 

2.80 

0.18 

1.00 

1.23 

54°Xe66 

1.16(10) 

2.19 

1.17(19) 

2.98 

0.96(17) 

3.82 

0.91(16) 

4.40 

2.70 

0.18 

0.91 

1.17 

^fXe68 

1.46(11) 

1.95 

1.41(9) 

2.44 

1.03(8) 

3.01 

1.54(10) 

3.33 

1.90 

0.13 

0.47 

0.86 

ifXeyo 

1.17(4) 

2.10 

1.52(14) 

2.76 

1.14(36) 

3.48 

0.36(5) 

3.94 

2.40 

0.16 

0.72 

1.04 

ifXer2 

2.20 

2.99 

3.83 

4.42 

2.71 

0.18 

0.93 

1.18 

^fXer4 

1.47(15) 

2.21 

1.94(20) 

3.03 

2.39(30) 

3.89 

4.50 

2.76 

0.18 

0.97 

1.20 

54°Xe76 

2.28 

3.23 

4.22 

4.97 

2.99 

0.19 

1.18 

1.34 

i"Ce74 

1.11(26) 

2.00 

1.51(87) 

2.53 

0.73(16) 

3.15 

0.47(12) 

3.51 

2.06 

0.14 

0.54 

0.91 

i"Ce76 

0.75(17) 

1.91 

0.27(10) 

2.35 

0.47(4) 

2.88 

0.07(2) 

3.17 

1.75 

0.12 

0.41 

0.82 


point have been observed within their chains. 54 ®Xe 74 was also proposed in Ref. na as a 
good candidate for the E(5) symmetry, while in Refs. [26l[^ where a sextic potential is 
used to describe the shape phase transition from spherical vibrator to y-unstable nuclei, 
gl'^Xero and s^^Xegs are indicated as being the most closest to the corresponding critical 
point. All the above examples suggest that the shape phase transition predicted by 
X(3)-Sextic has a strong experimental support. 
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Table 7. The same as in Table [6l but for the available experimental data [54117(1] 
of the nuclei 150,152gj^^ 152,154 qj^ 154,156i3y^ 172Qg^ 180-196p^^ 190jjg 

222Ra. 


Nucleus 

^0,4,0,2 

Po,6,0,4 

Po,8,0,6 

Pq,10,0,8 

Pi,0,0,2 

Pi,2,0,2 

Pi,2,0,4 

Pi,2,1,0 


1.47(39) 

2.30 

3.31 

4.35 

5.17 

3.06 

0.20 

1.27 

1.40 

If^dss 

1.62(9) 

2.22 

1.76(14) 

3.06 

1.69(30) 

3.94 

4.57 

0.54(4) 

2.80 

0.25(3) 

0.18 

0.28(14) 

1.00 

1.23 

^e°Nd9o 

1.56(4) 

1.66 

1.78(9) 

1.93 

1 .86(20) 

2.27 

1.73(10) 

2.43 

0.37(2) 

0.95 

0.09(3) 

0.08 

0.16(6) 

0.16 

1.38(112) 

0.64 

62°Sm88 

1.93(30) 

2.29 

2.63(88) 

3.28 

2.98(158) 

4.30 

5.09 

0.93(9) 

3.03 

0.20 

1.23 

1 nQ+O-70 
^■'7 'J_o.53 

1.37 

ei^Smgo 

1.44(2) 

1.75 

1.66(3) 

2.07 

2.02(4) 
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Before going to conclusions, some final remarks are worth to be made here in order 
to frame the present solution between other Bohr-Hamiltonian approaches. X(3)-Sextic 







Figure 4. The plot of the free parameter a as a function of the neutron number N 
for the isotopes of Ru (a), Xe (b), Nd (c) and Pt (d). 


is adding to the special class of rigid solutions, which implies that one or more variables 
are frozen. Most of these solutions have the advantage of being exactly determined, 
obtained in an analytical form, not depending on free parameters and giving good 
agreement with the corresponding experimental data. Of course there are also some 
drawbacks, the major one being that by freezing some degrees of freedom a part of 
information is definitively lost. A good example in that sense is the impossibility 
to describe the 7 band. An exception is the Z(4) type solutions [mUS], where the 
description of 7 band is still possible. Nevertheless, this type of solutions should not be 










Sextic potential for 'y-rigid prolate nuclei 


18 


seen as stand alone descriptions but rather as complementary approaches. As has been 
discussed above, these solutions can be useful in the understanding of the critical point 
phenomena, in the clarifying of the nature of some ambiguous states and of course in 
the generating of new development perspectives for the present held or other branches 
of physics, respectively. 

4. Conclusions 

A new solution of the Bohr-Mottelson Hamiltonian HH with a sextic potential ESI 
is proposed, this time for 7 -rigid prolate nuclei. The model is called X(3)-Sextic. The 
separation of variables is exactly achieved. The angular equation corresponds to the 
spherical harmonic functions, while that for the (3 variable with a sextic potential is 
reduced to a quasi-exactly solvable form [20]. Finally, the energy spectrum and the 
wave function are given in analytical form and moreover, up to a scale factor, depend 
on a single free parameter. In particular cases of the sextic potential, when (3“^ or (3^ 
cancels, parameter free solutions are obtained. As the properties of the states 0 ^ and 
2^^ are exactly determined, the model becomes useful in deciding whether or not some 
ambiguous states are the heads of the hrst two (3 bands. The comparison of the X(3)- 
Sextic results with those of previous solutions as X(3)-/3^ [IT], X(3)-/d"^ [18], X(3)-/5® [18] 
and X(3) [15] reveals a good agreement between them, establishing in this way that the 
approximations involved to describe states characterized by L > 4 are rigorous enough. 

Numerical applications are done for 39 nuclei presenting experimental data at least 
for the ground band and hrst (3 band, namely for 

146-150Nd, 150.152g^^ 152,154 q^^ 154,156i3y^ 1720s^ 180-196pt^ igOpg 222j^g^_ 

nuclei where the band head states of the hrst two (3 bands were uncertain, a selection 
has been done by choosing those of them which were closer to the predicted ones. From 
all these applications is evidenced not only the accuracy of X(3)-Sextic to reproduce the 
existing experimental data, but also its power of prediction. 

X(3)-Sextic describes very well only the (3 bands of some X(5) candidate nuclei, as 
^®°Nd, ^®^Sm, ^^'^Gd and ^^®Dy, acting as a complement of X(5) which describes very 
well the ground state band, the 7 bands and the head of the hrst [3 band. Therefore, a 
combination of the two solutions as in Ref. [37] by involving a rigidity control parameter 
seems to be a proper method for these nuclei. 

By studying the dependence of the energies on the free parameter, a hrst order 
shape phase transition is pointed out covering a region between 7 -rigid prolate harmonic 
vibrator and 7 -rigid prolate anharmonic vibrator. This behavior is also experimentally 
evidenced within the isotopic chains of Ru, Xe and Nd, where the best candidates for 
the critical point are ^°^Ru, and ^^®Nd, respectively. Other possible candidates 

can be also considered ^®®Pt, ^^^Os or ^^®Xe. 
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The eigenfunctions and the eigenvalues of the quasi-exactly solvable 

equation fl2.19p are given here in an explicit form for M = 0,1, 2. 

The case M = 0: 

dhy')= K°l =0- 

The case M = 1: 


(3.27) 


paid-)= 


1 - 


^n,L 2 

:y 


A, 


( 1 ) 

0,L 


The case M = 2 
A2) /„2 


= 2 
= 2 


‘2(yL + 3) 

Oi — \/ + 2{L + 3) 

a + \lo? 2[L + 3) 


, n e {0,1}, 


(3.28) 


pyUv^) = <1 

= 4a 


),{2) 


2ASi 


2(L + 3) (i + 3)(Ara - 8a) 


2 / 


, n E {0,1,2}, 


),( 2 ) 

\,L 


2i{-i + \/3)(8 + a^ + 2L) ^ 2i{i + V3)D{L, a) 


3i 


X(2) 

\,L 


= 4a + 


2)3D{L, a) 

2i(i + x/3)(8 + a2 + 2L) 2i{-i + ^/2)D{L,a) 


2)3D{L, a 


3i 


, (3.29) 


A 


( 2 ) 

2,L 


, 4(8 + a2 + 2L) , AD{L,a) 
4a ~\~ 1 ~\~ 


33D{L,a) 


3i 


D{L,a) = 9a + V3^/27a^ - (8 + a^ + 2L)3 ' , i = ^/^. 

Here are constants which are determined from the conditions that the ansatz 
functions fl2.18p to be normalized to unity with the integration measure dy. These 
constants can be determined either numerically or analytically by evaluating integrals 
of the form 






A + 1 1 a^ 


2 2 


(3.30) 


where f/ is a confluent hypergeometric function CU- It is necessary to point out that all 
solutions are real dehned due to the special relations realized between the parameters 
of the quasi-exactly solvable sextic potential. 
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